1 



Experimental review of unpolarised nucleon structure functions 

in ' 

O Vladimir Chekelian (Shekelyan) E 



X 



a Max-Planck-Institut fiir Physik (Werner-Heisenberg-Institut) 
Fdhringer Ring 6, 80805 Miinchen, Germany 
pL( ■ E-mail: shekeln@mppmu.mpg.de 
CO 

Recent results are reviewed on unpolarised structure functions from fixed target exper- 

^ ' iments at JLAB, NuTeV and from the HERA ep collider experiments HI and ZEUS. 
OO ' 

o 

O 1- Introduction 

CN 

The study of the structure of hadrons has been a powerful means for establishing and 
testing of the theory of strong interactions and for the determination of the partonic 
content of the nucleon. In 1969 the SLAC-MIT collaboration [ observed the scaling 
behaviour of the proton structure function in deep inelastic electron-proton scattering 
Q_i| (DIS). This observation established the quark-parton model (QPM) as a valid framework 
for the interpretation of DIS data and that different partonic processes can be expressed 
in terms of universal parton densities. The smallness of R = ctl/vt, the ratio of the 
cross sections from longitudinally and transversally polarised virtual photon scattering 
measured in DIS [Ej, provided the first evidence of the spin- 1/2 nature of the partons. 
In 1974 in fj,N interactions [ 4j and in subsequent neutrino and muon-nucleon scatter- 
ing experiments, the observation of scaling violation and the identification of partons as 
quarks and gluons has confirmed the field theory of quarks and gluons and their strong 
interactions, Quantum Chromodynamics (QCD). 

The structure functions in QCD are defined as a convolution of the universal par- 
ton momentum distributions inside the proton and coefficient functions, which contain 
information about the exchanged boson-parton interaction. At sufficiently large four- 
momentum transfer squared, Q 2 , when the strong coupling a s is small, a perturbative 
technique is applicable for QCD calculations of the coefficient and splitting functions. 
The latter represent the probability of a parton to emit another parton. Expressions for 
structure functions are determined by convolution integrals of appropriate sums over the 
densities of quarks of different flavours and gluons, which predict a logarithmic Q 2 depen- 
dence (evolution) of the structure functions. The perturbative QCD (pQCD) calculations 
are well established within the DGLAP formalism to next-to-leading (NLO) order in 
the strong coupling and have recently been extended to next-to-next-to-leading (NNLO) 
order 

In 2004 D. Gross, D. Politzer and F. Wilczek have been awarded the Nobel Prize for 
the discovery of asymptotic freedom. F. Wilczek in his brief commentary [|7] on the QCD 
foundational papers wrote that the "most dramatic" experimental consequence "regarding 
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the pointwise evolution of structure functions", namely "that the proton viewed at ever 
higher resolution would appear more and more as field energy (soft glue), was only clearly 
verified at HERA twenty years later". Indeed, the first measurements by the HERA 
electron proton collider experiments HI and ZEUS in 1992 [IHUHl revealed a steep rise of 
the proton structure function towards small x, the fraction of proton momentum carried 
by the parton. The present state of the art of these measurements [ [TUl HH E2J EH E] is 
shown in Figured] as a function of x at Q 2 = 15 GeV 2 (left figure) and as a function of 
Q 2 at different x values (right figure). 
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Figure 1. Proton structure function F 2 (x,Q 2 ) measured at HERA and in fixed target 
experiments as a function of x at Q 2 = 15 GeV 2 (left) and as a function of Q 2 for different 
x values (right). Curves represent different NLO QCD fits. 



Partons, i.e. quarks and gluons, enter differently into different structure functions. In 
deep inelastic scattering the cross section of the neutral or charged current (NC or CC) 
process can be expressed in terms of three structure functions. For example for neutral 
current process mediated by 7 or Z°-boson exchange, the expression is 

dxdQ2 = ^ [Y + F 2 (x, Q 2 ) - y 2 F L (x, Q 2 ) T Y_xF 3 (x, Q 2 )] , (1) 

where y = Q 2 /xs is the inelasticity, s is the center of mass energy squared, a is the fine 
structure constant and Y± = 1±(1 — y) 2 . The dominant contribution to the cross section 
is due to the proton structure function F 2 , which, in the framework of QPM, is related 
to the sum of the proton momentum fractions carried by the quarks and antiquarks in 
the proton weighted by the quark charges squared. The longitudinal structure function 
Fl, vanishing in QPM, is directly sensitive to the gluon momentum distribution in the 
proton. The structure function xF s depends on the valence quarks and is sizable only at 
large Q 2 , when Q 2 is comparable with the Z°-boson mass squared. 
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In this paper the recent JLab results on the transverse and longitudinal structure func- 
tions in the nucleon resonance region, the NuTeV measurements of the isoscalar structure 
functions F 2 (x,Q 2 ) and xF 3 (x,Q 2 ) in vFe interactions and the HERA results on struc- 
ture functions are reviewed as presented at the BARYONS 2004 conference. The HERA 
results include analyses of the inclusive NC and CC e p cross section measurements in the 
framework of perturbative NLO QCD, determination of the quark and gluon distributions 
inside the proton, results on the strong coupling obtained from inclusive and jet data and 
also measurements of the charm and beauty contributions to F 2 . 

2. Recent results from fixed target experiments 

The Jefferson Lab experiment E94-110 measured inclusive scattering of unpolarised 
electrons off a hydrogen target in the nucleon resonance region, in the range 0.2 < Q 2 < 
5.5 GeV 2 The longitudinal-transverse separation allowed the proton structure func- 

tions F 2 = (2xFi + Fl)/(1 + 4M 2 x 2 /Q 2 ), F\ (purely transverse) and Fl (purely longitu- 
dinal) to be extracted independently. Here, M is the proton mass. The measurements of 
2xFi and Fl are shown in Figure El as a function of x for different Q 2 values. The longitu- 
dinal component is found to be significant, both structure functions show resonant mass 
structures which oscillate around the lines corresponding to QCD fits. At Q 2 > 1.5 GeV 2 
the average x dependence of the resonance region is well described. The observed scaling 
relationship between resonance electroproduction and deep inelastic scattering, termed 
quark-hadron (Bloom-Gilman) duality [El], suggests a common origin for both kinematic 
regimes and brings additional information on the transition from soft to hard QCD. These 
data represent the first observation of duality in the separated transverse and longitudinal 
structure functions. 
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Figure 2. The transverse (2xF\) and longitudinal (Fl) structure functions from the E94- 
110 experiment (triangles) compared with the SLAC measurements (diamonds) and the 
NNLO QCD fits from Alekhin [Ej (dashed line) and MRST [EE] (solid line) with, and 
MRST (dotted line) without target mass effects included. The prominent resonance mass 
regions are indicated by arrows. 
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Figure 3. The NuTeV measurements of F 2 (x,Q 2 ) and xF 3 (x,Q 2 ) compared with the 
data from CCFR and CDHSW. The curve represents the fit to the NuTeV data. 



The NuTeV experiment has obtained high statistics samples of neutrino and antineu- 
trino CC interactions using the high energy sign-selected neutrino beam at Fermilab. 
Isoscalar neutrino-iron structure functions F 2 (x,Q 2 ) and xF 3 (x, Q 2 ) [E] determined by 
fitting the y dependence of the sum and the difference of v and v cross sections are shown in 
Figure El The measurements benefit from the precise determination of muon and hadron 
energy scales (0.7% and 0.43%, respectively) obtained using continuous calibration beam 
running concurrently with the data taking. At intermediate x (0.015 < x < 0.5) both 
F 2 (x,Q 2 ) and xF 3 (x, Q 2 ) are in good agreement with the CCFR results which are also 
shown in Figure El At lowest x (x = 0.015) and in the high x region (x > 0.5) the 
NuTeV data are systematically above the CCFR results by ~ 3 % and by up to ~ 20 %, 
respectively. The new measurements have improved systematic precision, and they have 
expanded into the previously unaccessible kinematic range of high inelasticity y. 

3. HERA results 

The kinematic coverage of the HERA structure function measurements in the x-Q 2 
plane is shown in Figure 0] (left). HERA extends the phase space of the previous fixed 
target measurements (also shown in the figure) by 2 orders of magnitude both in x and 
Q 2 . The full HERA range in x is essential for predictions for the LHC collider. 

At the HERA collider, the incident electron (or positron) energy is 27.6 GeV and the 
proton energy is 920 GeV (820 GeV till 1997), which correspond to a center of mass 
energy of 319 (301) GeV. In the first phase of data taking, HERA I, the HI and ZEUS 
experiments each collected about 100 pb _1 of integrated luminosity in the e + p mode and 
about 15 pb _1 in the e~p mode. The HI and ZEUS measurements of the single differential 
NC and CC e ± p cross sections dcr/dQ 2 [ H31 Q3J 1201 12H 123 are summarised in Figure H 
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Figure 4. Kinematic reach at the HERA collider and in the fixed target experiments 
together with the regions of jet production in pp scattering (left). The HI and ZEUS 
measurements of the NC and CC e^p cross sections as a function of Q 2 (right). 



(right). At low Q 2 ~ 100 GeV 2 the cross section of the CC process mediated by the 
W boson, is smaller by 3 orders of magnitude than that of the NC process, due to the 
different propagator terms. At high Q 2 « Mf , Myy the cross section measurements are 
approaching each other demonstrating the unification of the weak and electromagnetic 
forces. From a comparison of the NC measurements at highest Q 2 with the Standard 
Model expectation, a limit on the quark radius of ~ 10~ 18 m is obtained, proving a 
pointlike behaviour of quarks down to about 1/1000 of the proton radius. 

3.1. Partonic structure of the proton 

The measurements of the full set of NC and CC double differential e ± p cross sections at 
HERA allow comprehensive QCD analyses to determine the quark and gluon distributions 
inside the proton and the strong coupling constant a s (M§). 

Double differential NC e ± p cross section measurements at high x are shown in Fig- 
ure El (left). At low Q 2 the cross sections of the e + p and e~p interactions are essentially 
indistinguishable. At high Q 2 they depart from each other due to the different sign of 
the xF% contribution to the cross section (see eq. QJ). The HERA experiments determined 
the structure function xF 3 [ [T3J 122] by the difference of e + p and e~p cross sections. It is 
dominated by the 7 — Z interference term and depends on the valence quark density only. 
Double differential CC e^p cross sections are shown in Figure 03 (right) as a function of 
x for different Q 2 . They are sensitive to individual flavours in the proton, the e + p data 
to d and the e~p data to u quarks. At high x, the contribution from the valence quarks 
dominates the cross section and allows a local extraction of u and d densities. 

Using inclusive measurements of the NC and CC e ± p cross sections, HI and ZEUS 
performed NLO QCD fits [E3 |2I], which lead to a quark flavour decomposition of parton 
densities. In the ZEUS fit their jet data are used as well. The quark and gluon distribu- 
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Figure 5. Double differential NC (left) and CC (right) e ± p cross sections at HERA 
compared with the Standard Model predictions based on the CTEQ parton distributions [ 
. For CC, the contributions from u and d quarks are also shown. 



tions obtained in the HI and ZEUS fits are shown in Figure IBl (left). The results agree 
within the quoted error bands. They agree also with the parton densities from global 
QCD fits I23 which include not only the HERA but also the fixed target DIS data 
as well as data from other processes sensitive to parton distributions, such as Drell-Yan, 
inclusive jet production and W-lepton asymmetry in pp collisions. 

3.2. Low-i regime 

The steep rise of the proton structure function at low x indicates a new partonic regime 
governed by a rising gluon density. This high parton density regime might reveal novel 
QCD effects, e.g. gluon-gluon recombination. The QCD analysis of inclusive measure- 
ments leads to a gluon distribution which rises towards small x at Q 2 above a few GeV 2 , 
Figure |Hl (right). At the lowest Q 2 the rise of F 2 persists, while the gluon density flattens 
out at Q 2 ~ 2.5 GeV 2 and even becomes dangerously close to zero at Q 2 = 1 GeV 2 . This 
is not necessarily a problem in itself, since the gluon density from scaling violations is 
not an observable. A negative gluon density, however, would result in a distinct unphysi- 
cal prediction for Fl- A closer look into the low x region is presented below in terms of 
derivatives of F 2 in In x and In Q 2 and measurements of the longitudinal structure function 
/•/, 

In the double asymptotic limit, the DGLAP evolution equation can be solved and F 2 is 
expected to rise approximately as a power of x towards low x. A power behaviour is also 
predicted in BFKL theory . A damping of this rise would indicate the presence of novel 
QCD effects. A relevant observable for the investigation of the dynamics of this growth 
is the partial derivative of F 2 w.r.t. lnx at fixed Q 2 , A = — (<9 In 1*2(0;, Q 2 )/d\nx)Q2 . 
The high precision of the present F 2 data allowed HI to measure this observable locally [ 
l2*7j . The measurements are consistent with no dependence of A on a; for x < 0.01. 
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Figure 6. The valence, sea and gluon distributions as obtained from the HI and ZEUS 
NLO QCD fits to NC, CC and jet data (latter in ZEUS fit only) at Q 2 = 10 GeV 2 
as a function of x (left). The low x region is dominated by the gluon and sea quark 
distributions divided on the plot by a factor of 20. The gluon distribution from the ZEUS 
fit at Q 2 = 1, 5, 20 and 200 GeV 2 (right). 



Thus, the monotonic rise of F 2 persists down to the lowest x measured at HERA, and 
no evidence for a change of this behaviour such as a damping of the growth is found. 
The observed independence of the local derivatives in hi 2 at fixed Q 2 suggests that F 2 
can be parameterised in a very simple form F 2 = c(Q 2 )x~ x ^ * . The results for X(Q 2 ) 
obtained by HI and ZEUS are shown in Figure [7| (left). The coefficient c(Q 2 ) ~ 0.18 and 
the parameterisation X(Q 2 ) = a-ln(Q 2 /A 2 ) for Q 2 > 2 GeV 2 are consistent with pQCD 
analyses. At Q 2 < 1 GeV 2 the behaviour is changing, and, in the photoproduction limit 
(Q 2 ~ 0), A is approaching 0.08, which is expected from the energy dependence of soft 
hadronic interactions a to t ~ s a p(°)- 1 ~ s - 08 . 

Another important quantity in view of possible non-linear gluon interaction effects is 
the derivative (dF 2 /dln Q 2 ) x which is a direct measure of scaling violations. Its behaviour 
in a; is a reflection of the gluon density dynamics in the associated kinematic range. The 
derivative measurements are shown in Figure [7| (right) as a function of x for different Q 2 . 
They show a continuous growth towards low x without an indication of a change in the 
dynamics. The derivatives are well described by the pQCD calculations for Q 2 > 3 GeV 2 . 

Non-zero values of the structure function Fl appear in pQCD due to gluon radiation. 
Therefore, Fl is a most appropriate quantity to test QCD to NLO and especially to 
examine pathological effects related to a possibly negative gluon distribution. According 
to eq. [TJ the Fl contribution to the inclusive cross section is significant only at high y. The 
conventional way to measure Fl is to explore the y dependence of the cross section at given 
x and Q 2 by changing the center of mass energy of the interaction. Such measurements are 
not yet performed at HERA. The HI collaboration nevertheless could determine Fl from 
measurements at high y, i.e. small scattered electron energies down to 3 GeV. Various 
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Figure 7. Fitted values of X(Q 2 ) (right). The derivative (dF 2 /d\riQ 2 ) x as a function of 
x for different Q 2 (right). 



methods are used which attribute the observed decrease of the cross section at high y 
to Fl according to eq. H A summary of the Fl measurements by HI [ |2H] is shown in 
Figure |Hl (left). The results are significantly above zero everywhere, including the lowest 
Q 2 . They are compared with pQCD calculations and different phenomenological models 
showing that already at the present level of precision the measurements can discriminate 
between different predictions. Direct measurements of Fl at HERA can be performed 
only by reducing the beam energy and employing the highest y domain. 

3.3. Strong coupling a s 

Further insight into the structure of the proton is obtained in semi-inclusive processes 
such as jet or heavy flavour production which are calculable in pQCD. The jet data 
have been included into the QCD analysis by ZEUS [ EI] • They help to reduce the 
uncertainty of the gluon distribution at medium x around 0.1. In this fit both, the gluon 
density shown in Figure |H1 (right) and the fundamental coupling constant a s (M|) were 
determined. The HI collaboration has also performed a QCD fit [ ITT] devoted especially 
to the determination of the gluon distribution and a s . A novel flavour decomposition 
of F 2 used in this fit allows to reduce the number of parton distributions to just two 
combinations (apart from the gluon density): "valence-like" and "sea-like". The number 
of data sets can thus be reduced to a minimum, the HI cross section data covering lowx 
and the BCDMS fip data covering the region of large x. Without HERA data the gluon 
distribution from the fit to the BCDMS data becomes unrealistically flat at low x and 
drives a s (M|) to about 0.110. The HERA data pin down the gluon density, and, as a 
result, a s (M|) moves to 0.115 if both data sets are used. 

The most precise HERA measurements of the strong coupling constant a s (M|) in 
inclusive DIS are shown in Figure |H] (right) together with various a s (M§) determinations 
from jet production at HERA [I2H]- The average a s (M|) value in NLO from HI and 
ZEUS is 



a s {M 2 z ) = 0.1186 ± O.OOll(exp) ± 0.005(theory). 



(2) 
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Figure 8. Summary of the Fl measurements by HI at fixed photon-proton center of mass 
energy W = 276 GeV, W ~ ^/sy (left). Measurements of the strong coupling constant 
in inclusive DIS and jet production at HERA (right). 
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The experimental error is small, about 1%. A limitation of the precision of a s (M|) 
arises from the associated theoretical error calculated using the ad hoc convention of 
varying the renormalisation and factorisation scales by a factor of 2. In forthcoming 
exact NNLO analyses the scale dependence will be considerably reduced. With theoretical 
and experimental progress in DIS, it is expected to pin down the least well- measured 
fundamental constant a s (M|) to an accuracy of better than 0.001. 

3.4. Heavy flavour production 

The dominant heavy flavour production process at low x at HERA is photon-gluon 
fusion, in which a pair of charm or beauty quarks and antiquarks is produced. The mass 
of the heavy quark provides a hard scale and makes it possible to apply pQCD techniques 
for cross section calculations even in photoproduction (Q 2 w 0). The interplay between 
different hard scales, mass and pt of the heavy quark and Q 2 , is one of the questions which 
heavy flavour studies need to resolve. 

Open charm production at HERA has mainly been studied by reconstructing D* mesons 
which originate from the hadronisation of charm quarks and decay into D°tt — > K7nr. 
The measured contribution of charm to the proton structure function [EDI EH] is shown in 
Figured (left). The charm contribution is large and increases from w 10% at Q 2 = 2 GeV 2 
to ~ 30% at Q 2 = 500 GeV 2 . Silicon detectors surrounding the interaction region provide 
lifetime tagging for heavy flavour physics at HERA. Using this powerful tool the HI 
collaboration determined the cross sections for charm and beauty production [|32] shown 
in Figure El (right). The vertex tagged HI charm data are consistent with the ZEUS 
D* data and complement the data obtained at lower Q 2 . The beauty contribution is 
about 1/10 of the charm contribution and amounts to about 2% of F2. This is the first 
measurement of the beauty structure function The data are well described by pQCD 
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Figure 9. Fractional contributions of charm and beauty quarks to the proton structure 
function F 2 . The curves correspond to predictions within the HI and ZEUS QCD analyses. 



calculations both for F 2 CC and F$ . The accurate measurement of these structure functions 
is important for the forthcoming LHC data, because their contribution is expected to be 
much increased for the scale relevant for the LHC. 

3.5. First results from HERA II 

Data taking of the second, high luminosity phase of the HERA program, HERA II, 
started in October 2003. After the upgrade of the collider the specific luminosity is 
increased by about a factor of 3. This was achieved by placing strong super-conducting 
focusing magnets inside the HI and ZEUS detectors, close to the interaction point. In 
2003-2004 HERA has delivered about 100 pb" 1 . A major success at HERA II is the 
longitudinal polarisation of the positron beam (P), typically about 40%. This new feature 
was exploited by HI and ZEUS to measure polarised NC and CC e + p cross sections [ 
EHl Ell- The measured total CC cross sections are shown in Figure ITU1 together with 
unpolarised HERA I data (P = 0). The measurements confirm the predicted linear 
dependence of the CC cross section on the longitudinal polarisation of the positron beam. 
The e + p cross section is expected to vanish at P = — 1 unless right handed weak currents 
exist. The result of the combined HI and ZEUS cross section measurements, linearly 
extrapolated to P = -1, is 0.2 ± l.8(stat) ± l.Q(syst) pb (x 2 /dof = 5.4/4), which is 
consistent with absence of right handed weak currents. 

4. Outlook 

For more than three decades unpolarised structure function measurements in DIS are 
providing crucial experimental input to establish the quark parton model, to develop 
QCD and to determine universal parton distribution functions. The area is still very 
active and provides many results from fixed target experiments and, in the past decade, 
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Figure 10. Measurements of the total CC e + p cross section as a function of the positron 
beam polarisation P. The published HERA I results correspond to P = 0. The first 
preliminary HERA II measurements are made at P = —0.40 and P = 0.33. The lines 
indicate the Standard Model expectation and a linear fit to the data. 

from the unique ep collider HERA. The second phase of the HERA program with increased 
luminosity and upgraded detectors just started. The rich physics program of HERA II 
has the goal to collect about 1 fb _1 of data, including a special run with lower proton 
beam energies for direct measurements of Fl at low x. The end of HERA, foreseen for 
mid 2007, will definitely not allow to complete the program and will enforce compromises 
on the physics outcome from HERA. The field has a large potential in a long term future. 
Extensions of the present HERA program to study neutron structure in ed collisions and 
for measurements at low x and Q 2 with improved precision have been proposed recently [ 
13*5] . The HERA (or Fermilab) proton ring could serve for ep collisions with the future 
International Linear Collider at the new energy frontier [ I36j . Physics opportunities and 
the accelerator and detector options for a future Electron Ion Collider, EIC/eRHIC at 
BNL [EZj, are under discussion. 
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